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ABSTRACT: In the Villain approximation of D = 2 + 1 compact QED, 
the monopole part of the partition function factorizes from the Dirac string 
part, which generates the photon propagator. Numerical experiments in exact 
compact QED confirm this result: photon mass pole My, originally nonzero, 
is insensitive to monopole prohibition but almost vanishes if Dirac strings are 
prohibited. 



I. Motivation 



Lattice calculations [1,2] have evaluated the mass poles Mg^^y^ of gluon, 
quark, photon, and electron propagators in LGT(lattice gauge theory). More 
recently, "photon" propagators in abelian-projected SU{2) LGT were stud- 
ied [3]. Remarkably the SU{2) photon mass, as in cQED3_|_i (compact QED in 
3 + 1 dimensions) [2], is nonzero in the confine phase and vanishes in the decon- 
fine phase. Since computing these masses, which are gauge variant [4], requires 
gaugefixing to either Landau or maximal abelian gauges it is unclear whether 
they are indicators of real physics, lattice regularization artifacts, gaugefixing 
artifacts [5], or some combination thereof. 

The abelian-projection potentially provides a concrete framework for in- 
vestigating the relationship between LGT gluon masses and confinement via 
abelian-projected photon masses [6]. This possibility beckons us to clarify 
the relationship between nonperturbative photon masses and confinement in 
cQED. To this end, this Letter focuses on cQED2_|_i, a QCD-like model with 
local gauge invariance, chiral symmetry breaking, and area-law electron con- 
finement. cQED confines because quantum monopole fluctuations [7] restrict 
electric flux to Abrikosov tubes of width A in a dual Meissner effect [8]. Recent 
LGT simulations resolve a finite London penetration depth A which, following 
Landau- Ginzburg theory, can be identified as the inverse mass of an effective 
gauge potential Ajf . 

As described below, the nonperturbative zero-momentum photon propaga- 
tor mass pole My in cQED2_|_i is also resolvable and nonzero. The observation 
of two mass scales, A and My, raises the following question: How, if at all, is 
Ajf related to primary photon field and A to My! Are they equivalent? 
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Indeed, as shown below, A and My in CQED2+1 are not blood relatives. The 
gauge-invariant quantity A arises from quantum monopole fluctuations while 
My arises from gauge- variant Dirac strings. 

II. Monopoles, Gaugefixing, and Dirac Strings 

Let dphx = h x+ ^ - h x , d*h x = h x - h x _^, and □= d*d^. Link and 
plaquette angles in cQED, 

0^ = A^ + 2Trn^, -tt<A^<tt, n^Z (I) 

F^[A] = d^Ay - dyAp, (2) 

depend only on photon A^. The cQED action S c = /3^ /U < iy (l — cos F^) is 
invariant under local gauge transformations 

89 H = -d^u X} SA^ = (Ap - d^^Modf-Tr, tt) - A^. (3) 

While plaquette exp(iF^) is gauge-invariant, a gauge transformation inducing 
unequal shifts of on the four links of F^ v shifts F^ v by a 2tt multiple. F^ v 
decomposes into a gauge-invariant part £ [— 7r,7r) and a gauge- variant 

integral kink N^ v £ Z such that 

Ffj, v = (0 + 2vriV)^, (4) 

SF^ = 2ttSN^ = 8^(69 - 6A) V v). (5) 

Because cos F^ v = cos ©^j, in S c , required since N^i> is gauge- variant, ^F 
would be a misleading approximation to S c - More valid is Villain's periodic 
gaussian approximation S c —> where 

e SY ^^ e -fE,™-2^) 2 , (6) 
{N} 
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E{#} = n^E^^-oc^-^V^W)' and p2 = Hn,v F lv The sum over 
{N} in (6) is necessary to maintain gauge invariance under (5). Within an N^ v 

potential well the Maxwell equation is d^F^ v = K v = 2n '<9^iV ^ v . Conserved 

current acts as a virtual source of electromagnetic charge. Quantum N^ v 

fluctuations can lead, via the Maxwell equation, to a nonzero effective My if 

(O^n = M* (0^) N + ■■■ (7) 

In (7) ( )jy refers to the Boltzmann average over N^ v and operator O ^ serves 
to absorb D = 2 + I rotational invariance. 

While v is gauge variant, spatial combinations of them form gauge- 
invariant structures which influence as follows. According to the Hodge- 
DeRham theorem 

®nv = e ^ad*<p + dfj,a v - d^a^, (8) 

N^ v = e^ va d%m + d^l v - d v l^ (9) 

where </>, G ( — 00,00) and m, G Z. <f> and m are invariant under (3), 
transforms like 0^, and like (A — a)^l2ir. (and similarly N^ v ) 

has 3 independent polarizations while <f> and are 4 functions because 
is invariant under Sa^ = —d^u x - In vector notation Eq. (4) becomes B = 
H + 2irff where the total B } physical H } and Dirac string ff magnetic(actually 
electromagnetic) fields are 

BeVxA, H = V(f) + Vxa, rf = Vm + V x K (10) 

It will be advantageous to recast ff in terms of its divergence and curl 

q = V • if = Dm, p = V X ff = V(V • I) - Ul. (II) 

Since V • B = by (10), V • H = —2irq, that is, q causes dislocations in 
the physical field H . By tautology q is the magnetic monopole density, gauge 
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Figure IB. The same lattice after Landau gaugefixing. While monopole 
positions do not change, there are more Dirac strings. 
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invariant since m is gauge invariant. p } a continuous current wrapping around 
jf, is gauge- variant. 

In general, kinks occur either in monopoles, Dirac strings connecting a 
monopole antimonopole pair, or Dirac string loops. Loops can either be homo- 
logically trivial or toroidally wind around the periodic boundaries. Monopole 
charge density q is gauge-invariant but the number of string loops and the 
length and shape of all strings vary with gauge. Figure \A shows a monopole 
and Dirac string arrangement in a representative (3 = 1.0, periodic 5^ non- 
gaugefixed lattice. Figure IB shows the same lattice after 5000 Landau gauge- 
fixing sweeps. There are typically more Dirac strings after Landau gaugefixing 
than before, as depicted. 



III. Factorization of Photon Propagator from Monopoles 

Upon adopting the condition V • / = and ignoring Laplacian zero modes, 
Eqs. (4)-(ll) constitute 1-to-l variables changes {N} — > {m, /} — > {q } p} — > 
{(f), a} where, if nA X; y = -8 x ,y, d*A^ v = 0, and uA^ y = -8^ )V 8 X) y, then 
(f) = 27T ^x, y q y and = - 2ir Y, y ^yPvy- Following (4), (6), and (8) 

Z V C = [ e~ S c = Z m [0] x Z Al [0], (12) 
J A 

Z m \i\ = Y^e^ Cxqx - 2ir2p ^,y qxMx - y)q \ (13) 
{?} 

Z A l[0]= I e"fS^ 2 W 5 afi e (-00,00). (u) 
J a 

Eq. (12) relies on the quadratic character of the Villain approximation and 
e pv\^pv^\ c t > = Oj 0{®^) corrections from the cos 0s of S c would ruin fac- 
torization. (14) says a is a massless noncompact vector particle. 
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The dilute gas expansion and occupation number resummation over q £ 
{0,±1} of Z m in (13) yields [7] 

Z m [iW /e-A^^ $ -^ 2 - 2A " 2cOS<I) (15) 

where A^ = 27r^/3e _ ^ 7r ^A(0)_ Scalar $ is semiclassically identified with cf> in 
(8) because, since Y^x^qx = ^Ei^^^! comparing (13) to (15) yields 

i {d ^ )m = hwh-** = 2^ $) *" (16) 

Following (10) this identification yields 

(V-H) Alrn = -2*(q) m = 0, (17) 

(HyHx)Alm = ((V x 3)„(V x a) x ) Al + ^ J 2Z ™ . (18) 

If /3/A << 1, (17) and (18) are reproduced by an My = free photon A e ^ 
with H eS = V x A eff , that is, V • H eS = and (HfHf) eS « {H y H x ) Aim- 
The latter relation relies on the masslessness of a in (18), a consequence of 
(14), and cos $ — > — $^/2 in (15). is the massive Landau- Ginzburg photon 
and A the London penetration depth. 

The A propagator is generated by Z^y[J], defined by adding J ■ A to the 
action in (6), which does not affect factorization (12). Thus the A mass My 
has nothing to do with monopoles q and, hence, My is not directly related 
to A. A remaining logical possibility, which we do not exclude, is that A and 
My are coincidentally related via the statistical mechanics of monopoles and 
p loop gasses. A, unlike a, may be massive because J ■ A breaks the pure 
a-dependent form of Z^y[0]. Manipulations like those leading to (13) yield 

Z Al[J] = f ^eS^^+^^-f^-^E^-A-^ 
A {p} 
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Figure 2. (A)cQED with action S c ; 

(B) cQED with monopoles forbidden; 

(C) cQED with kinks forbidden; 

(D) quadratic approximation S c — > S^, kink-creating gauge transformations allowed; 

(E) quadratic approximation S c — > S^, kink-creating gauge transformations forbidden. 
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Zai is the partition function of a Coulombic p loop gas, "loop" since V • p = 
0. Interestingly p is a mixed state in the gas since for My to be nonzero 
Y2y yi A^Aq^, (py^a Py'fi)p must have a negative norm massless mode to 
cancel the a pole and an independent My mode. Note that even though 
photons are uncharged they suffer confinement since, heuristically, the "ad- 
joint" Wilson loop obeys i oop A l) ~ (ll/e loop sm ^) « Re (ILe loop e ^); 
where cross terms are suppressed by gauge invariance. 

IV. Numerical Experiments 

We have shown that the propagator decouples from monopoles in the 
Villain approximation and, accordingly, My is not a consequence of monopoles. 
Numerical experiments summarized in Figure 2 support this result in exact 
CQED2+1. (3 X My, a dimensionless number in D = 2 + 1, is the log of 
the ratio of successive p = photon propagator timeslices. The numer- 
ical photon operator and gauge condition are = sin#^ and d^S^ = 0. 
The computation is based on 500 configurations on 17 X 19 lattices. The 
first configuration — independently for each f3 value — is thermalized by 500 
forty-hit, 40%-acceptance Metropolis sweeps and 5000 checkerboard gaugefix- 
ing sweeps. Configurations thereafter are separated by 5 forty-hit Metropolis 
sweeps and 5000 checkerboard gaugefixing sweeps. Errors are jackknife sigmas 
based on 10 450-configuration subaverages. Configurations 1 — 50 are omitted 
from the first subaverage, 51 — 100 from the second, and so forth. 

"A" in Figure 2 refers to cQED in Landau gauge; "B" to Landau gauge 
cQED with monopoles prohibited [9]; "C" and "E" to cQED with kinks pro- 
hibited. Monopoles prohibition is implemented by starting with the 6^ = 
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configuration and linkwise forbidding monopole-creating updates with the in- 
sertion of n{x} ^qfi m t° ^ ne link measure. Kinks are prohibited either by 
inserting ri{#} ^ ne measure("C") or by replacing cos F^ v in S c 

with S E = f^("E") where (2) defines F^. 

Landau gaugefixing, which cannot change q } proceeds normally in B. Since 
maintaining kink prohibition disallows kink-creating gauge transformations 
which may be necessary to achieve Landau gauge, a good Landau gauge is not 
achieved in C and E. Nontheless the photon propagator is easily resolved in all 
cases and, in agreement with factorization (12), My is relatively insensitive to 
monopole prohibition but very sensitive to kink prohibition. The small residual 
My in C is attributable to the (9(0^) terms in S c which ruin factorization 
(12). "D" is also based on the action S^. Unlike in E, kink-creating gauge 
transformations are permitted during Landau gaugfixing in D. From the Sjt; 
standpoint D, corrupted by action-changing kink creation and annihilation, 
is gauge mequivalent to E. The difference between My in D and E, gauge 
equivalent from the S c viewpoint, is an indication of how much kinks generated 
by the Landau gaugefixing algorithm contribute to My. The smallness of My 
in D suggests that the kinks responsible for My in A are present in the pre- 
gaugefixing configurations and not created during Landau gaugefixing. 

My and kink-number density p } two gauge- variant quantities, are corre- 
lated. In cases A-E at f3 = 1.8 in Landau gauge, pj± = .41(.01), pfj = 0, 
PD = .23(.004), and ~ 10 - ^. Since the (3 = 1.8 monopole number den- 
sity is 8.0(1.1)10 - ^, forbidding monopoles doesn't change the kink density and 
PB = PA- I n t ne same f3 range the Axial gauge pj± is smaller than the Landau 
gauge pa by ~ 50% and My smaller by ~ 70%. 
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